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Gauge bosons associated with unbroken gauge symmetries, under which aU standard model fields 
are singlets, may interact with ordinary matter via higher-dimensional operators. A complete set of 
dimension-six operators involving a massless U{1) field, 7', and standard model fields is presented. 
The n 67' decay, primordial nucleosynthesis, star cooling and other phenomena set lower limits 
on the scale of chirality-flip operators in the 1-15 TeV range, if the operators have coefficients 
given by the corresponding Yukawa couplings. Simple renormalizable models induce 7' interactions 
with leptons or quarks at two loops, and may provide a cold dark matter candidate. 
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Experiments have established the existence of 12 gauge 
bosons. Ahhough gauge invariance requires the gauge 
bosons to be massless, three of these {W^,Z) are mas- 
sive due to spontaneous symmetry breaking, while eight 
other (gluons) are confined in massive hadrons. The only 
massless spin-1 particle discovered so far is the photon. 

This letter addresses the possibility that massless 
gauge bosons other than the photon exist. If any stan- 
dard model field is charged under a new unbroken gauge 
symmetry, the requirements of fermion mass generation 
and gauge anomaly cancellation force the new symme- 
try to be U{1) B-L, and the nonobservation of long range 
forces other than electromagnetism and gravity imposes 
an extremely severe constraint on the new gauge cou- 
pling, gs-L ^ m„/Mpi « 10~^^, where to„ is the neu- 
tron mass, and Mpi is the Planck scale (neutrino screen- 
ing cannot relax the bound 0). 

It is possible, however, that standard model fields are 
neutral with respect to a gauge group, but nevertheless 
interact with the new massless gauge bosons. This fol- 
lows from the existence of gauge invariant operators, of 
mass-dimension six or higher, which involve gauge field 
strengths. Such interactions have a natural decoupling 
limit: rather than forcing a dimensionlcss parameter to 
be extremely small, the experimental limits translate into 
a lower bound on the mass scale that suppresses the 
higher-dimensional operators. 

In the case of an Abelian gauge symmetry, U{l)p, un- 
der which all standard model fields are singlets, there is 
a single renormalizable term in the Lagrangian that in- 
volves both the U{l)p gauge boson and standard model 
fields: a kinetic mixing cqP^'^ B^^, where S'*" is the 
U{1)y (hypercharge) field strength, Pt'" is the U{l)p field 
strength, and cq is a dimensionlcss parameter. One might 
worry that this kinetic mixing would induce dimension- 
four couplings of the new gauge boson to all standard 
model fermions, but this is not the case "s]: the ki- 
netic terms of the two U{1) fields can be diagonalized 
and canonically normalized by an SL{2, R) transforma- 
tion, and their ensuing global 5*0(2) symmetry allows the 
identification of the linear combination of J7(l) fields that 
couples to hypercharge as the new hypercharge gauge 
boson. The orthogonal combination, referred to as the 



"paraphoton" in Ref. Q and denoted here by 7', does 
not have any renormalizable coupling to standard model 
fields. 

The leading interactions of the U{l)p field with stan- 
dard model fields are given by dimension-six operators, 
so they are suppressed by two powers of the mass scale 
M where the operators are generated. The terms in the 
Lagrangian of this type which involve fermions are given 
by 

+ ta'^'^CeHeR + h.c.) . (1) 

The notation is as follows: qljIl are quark and lepton 
doublets, URjdii are up- and down-type S'C/(2)-singlet 
quarks, en are electrically-charged S'C/(2)-singlet leptons, 
and H is the Higgs doublet. An index labeling the three 
fermion generations is implicit. The 3x3 matrices in 
fiavor space, C„,Cd,Ce, have complex elements which 
are dimensionlcss parameters. 

If "right-handed neutrinos" (gauge singlet fermions) 
are present, there is an additional dimension-six opera- 
tor involving P'"', as well as a lepton- number violating 
operator of dimension five: P^^N'^Ra^^ Nr. Otherwise, 
the leading 7' interactions with neutrinos are given by a 
dimension-seven operator, P^^Hl'^^a^^''' HIl. 

Eq. Q displays a complete set of dimension-six oper- 
ators involving U(l)p gauge fields and standard model 
fermions. For example, operators which involve the dual 
field strength, P^j, = Sfj^^xrP^'^ , can be reduced to Eq. ^ 
using the identities 

P^.a^'' = -2P^,(7^''75 , 
P^.^rO-' - P^, {a^'^Ip - 21^0'') 75 , (2) 

and field equations such as ilpqL ~ ^uHur + XdHdji, 
where \u,d are Yukawa coupling matrices, and the covari- 
ant derivatives D'^ contain only standard model gauge 
fields. Furthermore, using field equations and integration 
by parts, and ignoring operators of dimension higher than 
six, one can show that any chirality-preserving operator 
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can be written in terms of operators For example, 
the following identity is valid up to a total derivative: 

(3) 

In addition to the interactions with quarks and lep- 
tons shown in Eq. (Q, the U{l)p field has purely bosonic 
interactions described by dimension-six operators: 

H^H (ciB^, + hB^, + C2P^, + hh'^ P""" ■ (4) 



M2' 

These renormalize the U{1)y x U{l)p gauge couphngs, 
and include vertices with paraphotons, Higgs bosons, 
photons, and Z bosons. The dimensionless parame- 
ters ci,2, ci^2 are real. Other operators vanish {e.g., 
P^,^Df'H'^D''H and P^^G^G""'', where G'"' is the gluon 

field strength) or are total derivatives {e.g., P^i^G^G^'^). 

The dimensionless coefficients of the dimension-six op- 
erators can in principle have any value consistent with an 
effective theory description, namely below ~ 47r. How- 
ever, the operators flip chirality, and although the 
origin of flavor structure in the standard model is not 
yet known, the elements of G„,Grf and Ge are expected 
to be of the order of or smaller than the corresponding 
Yukawa couplings. 

Eq. 1^ is written in the weak eigenstate basis, where 
the standard model Yukawa couplings are flavor nondi- 
agonal. In the mass eigenstate basis, obtained by act- 
ing with different unitary 3x3 matrices on the left- 
and right-handed fields, the chirality-flipping operators 
change their flavor dependence: G/ ^ Gj = u[CfU{^\ 

where and [/^ are the unitary matrices that diago- 
nalize the masses of the / = e,u,d fermions. Note that 
Ul^'UI is the CKM matrix, wher eas the relation between 
the neutrino mixing angles and J7£ is looser. The inter- 
actions of mass-eigenstate fermions, /', with the U{l)p 
field appear in the Lagrangian as follows: 



Vh 
M2 



P^u [7^'^"^ (ReG} + z ImG} 75) /'] 



(5) 



These terms proportional to the vacuum expectation 
value of the Higgs doublet, Vh ~ 174 GeV, represent 
magnetic- and electric-like dipole moment operators. 

There are various phenomenological constraints on the 
7' interactions. First, the successful predictions of pri- 
mordial nucleosynthesis limit the number Ag^ of effec- 
tive relativistic degrees of freedom contributed by new 
particles that are in thermal equilibrium at a temper- 
ature Tbbn ~ 1 MeV. The maximum value of Ag^ is 
often expressed in terms of the maximum number of ad- 
ditional neutrino species allowed by the data: Ag'!^'^^ — 
(4/7) AiV™^''. The data on light element abundances ex- 
hibit some inconsistencies which translate into an uncer- 
tainty on AN^'^''. In Ref. 0, it is found that AN'J'"'' = 
0.6 at the 2a level, while in Ref. Q the 2a and 3a con- 
tours in the N^, versus baryon density plane extend to 
^^max Q_2 and 0.5, respectively. At any rate, the 



bound does not allow the two degrees of freedom of a 
paraphoton in thermal equilibrium. Thus, the parapho- 
ton must decouple at Ty> > Tbbn, so that during nucle- 
osynthesis it contributes an effective number of degrees 
of freedom |Q| 



A5,(Tbbn) = 2 



g* {T] 



BEN) 



-,4/3 



5*(Ty) 



(6) 



Given that g*(rBBN) = 43/4, the lower bound on the 
number of relativistic degrees of freedom at freeze out is 



g*{T^') > 27.5 X (AAr™'''^)-3/4 



(7) 



The number of degrees of freedom in the standard 
model is — 247/4 just above the temperature of the 
QCD phase transition, Tqcd = 150- 180 MeV 0, and 

= 69/4 just below Tqcd- Hence, Ty = 180 MeV is al- 
lowed as long as AN^J^^^ > 0.34. The freeze-out temper- 
ature would have to be higher by an order of magnitude, 
Ty > rur, a AiV™^'' « 0.25. Note that the bound on 
^jY^ax f].Qjjj a fit to the cosmic microwave background 
does not apply to the paraphoton, because at the time 
of recombination = 2 rendering the effective number 
of 7' degrees of freedom negligible [smaller by a factor of 
(43/8)4/3 than at Tbbn]. 

The interaction rate of the paraphoton with the ther- 
mal bath is given by 



r(T) = ^^^-3 



(8) 



where C(3) « 1.202, and {a^') is the thermally averaged 
cross section for 7' interactions with the standard model 
particles that are in thermal equilibrium. At tempera- 
tures just above Tqcd, the thermal bath includes light 
quarks {u,d,s), leptons (e,/x), photons, gluons and neu- 
trinos. The dominant 7' interactions involve the heaviest 
fermions, namely fi and s. A detailed study of the 7' 
decoupling would entail solving a set of coupled Boltz- 
mann equations, but for the purpose of estimating the 7' 
couplings in terms of T^' it suffices to make a simple ap- 
proximation: the 7' interaction rate, r(T^/), equals the 
expansion rate of the universe at freeze out. 



r(T-y 



rp2 

1' 



PI 



2^ 
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g*{T^') 



1/2 



(9) 



The 7' annihilation due to interactions with muons 
proceeds through the processes shown in Fig.l. The para- 
metric dependence of the annihilation cross section is 



{ay) 



M4 



(10) 



where = |(Gg)22| w/i/?7ip ^ 0(1) sets the strength of 
the muon-paraphoton interaction, and a is the fine struc- 
ture constant. Thus, the constraint on T^i results in a 
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FIG. 
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1: Paraphoton annihilation via muon pair production 



and Compton-like processes (7'/^ 



■ 7M 



limit on the effective mass scale of the muon-paraplioton 
interaction, 



M 



3.9 TeV X [g*{Ty)] 



7' 



1 GeV 



1/4 



(11) 



The Ty^ dependence dampens the sensitivity of the limit 
to the approximations used here. Using g*{T^') = 255/4, 
which includes the 7' degrees of freedom, and T^i ^ 180 
MeV, Eq. (mi gives M/^ Z 1.5 TeV. The Y^/ ss, 
7's 7s and 7's 7s processes at Ty impose a limit 
stronger by a factor of approximately {\/3ms/m^y^'^ « 
1.2 on M/^c7, where Cs = \{C'j)22\vh/ms. 

Another phenomenon that can be affected by parapho- 
tons is star cooling. One may use the studies of axions 
in order to derive the energy loss in stars due to 7' emis- 
sion 0. The axion is similar to the paraphoton: they are 
bosons with derivative couplings to fermions. The axion 
though has spin zero rather than one, so that the energy 
loss from stars is twice larger for 7' emission than for 
axion emission when the effective couplings are equal . 
The energy loss due to electron-paraphoton interactions 
is proportional to the square of ge-y' = Ace'in'^/M'^, where 
Ce = |(Cg)ii|w/i/me. The limit on 7' emission through 



Bremsstrahlung, such as e -I- He 



^He + 7', 



from the core of red giant stars (9| requires gl^,/{ATr) < 
2.5 X 10^27^ go ^j^g^^ M/y/c; ^ 3.2 TeV. Compton-like 
scattering, 76^ — > 7'e~, in horizontal-branch stars sets a 
slightly weaker limit, ^ 1.8 TeV. 

The neutrino signal from supernova 1987A also lim- 
its the cooling through 7' emission, setting a bound on 
the coupling of the paraphoton to nucleons. The leading 
nucleon-paraphoton interactions are similar to Eq. |SJl: 



M2 



(12) 



The CV form factors are of the order of (C^)ii or 
(C4)ii. For example, QCD sum rul es g ive ImCV ~ 
0.2Im(4C^ - C'Jn for the neutron [ij. The effec- 
tive nucleon-paraphoton coupling, gj\f.yi, is proportional 
to the nucleon mass: gj^^' « icnmdiTin/M'^, where 
c„ = \Cj^\vh/md ^ 0(1). Requiring that the super- 
nova was cooled predominantly by neutrinos implies 



gAT-r' 



^ 2 



10 



-10 



so that M/^ 7 TeV. There 



is also a range of larger nucleon-paraphoton couplings 
(2 X 10^^ ^ g^f-y' ^ 7 x 10^^) allowed by the supernova 
1987A signal: the paraphotons are there too strongly 
coupled to nucleons to escape easily the supernova, and 
too weakly coupled to produce a signal in the detectors 
that measured the neutrino signal. However, that range 
is not compatible with the nucleosynthesis constraint on 
quark-paraphoton couplings. 

The long-range forces induced by paraphoton exchange 
between chunks of ordinary matter are feeble. The inter- 
actions of nonrelativistic electrons or nucleons with the 
paraphoton are spin dependent. Given that the average 
spin of macroscopic objects is very small, the lower limits 
on M / y/c^ or M/ y/c^ imposed by measurements of long 
range forces are evidently weaker than the TeV scale. 

The magnetic- and electric-like dipole moments of the 
t quark could be probed at the Tevatron, LHC or a linear 
e+e" collider if the dimensionless coefficient |(Cu)33| is 
of order unity and the mass scale M is of order 1 TeV. 

The flavor off-diagonal interactions included in Eq. (O 
are constrained by various flavor changing neutral cur- 
rent processes. The most severe limit comes from the 



/i 67' decay, and depends on Ce^ 
comparison of the decay width, 

r (m ^ 67') 



\iCi)i,\vh/m,,. A 



(13) 



times the measured muon lifetime, 3.3 x 10^ eV~\ to 
the experimental limit shown in Fig. 2 of Ref. [lll| . 

Br (/i+ e+A) < 3 x 10^^ for any massless particle X, 



leads to 



M 



^ 15 TeV 



(14) 



Flavor-changing neutral current processes in the hadron 
sector are also affected by 7' emission. Note though that 
the TT~^j' decay is forbidden by angular momen- 

tum conservation, while the experimental constraints on 
other processes such as — s- pj' or B ^ K*y , where 
7' carries missing energy, are quite weak. 

In the presence of kinetic mixing of the U{l)p x U(1)y 
fields, the operators contribute to the magnetic 
and electric dipole moments of the quarks and leptons. 
Hence, there are constraints on the products of the ki- 
netic mixing parameter cq and the coefficients of opera- 
tors 1^. For example, the magnetic moment of the muon 
is shifted by ~ 4coc^(m^/M)2/e, and requiring that this 
is less than 10"^ gives M/^Jc^ Z 12 TeV. 

The bosonic operators l@| 
boson decays, h — s- 77', h 



lead to non-standard Higgs 
Zy and h — s- 7'7', with 7' 
behaving in detectors as missing transverse energy. If the 
mass scale M is not much higher than the electroweak 
scale and the parameters ci and ci are not much less than 
unity, a decay h 7^t could produce striking signals 
in collider experiments. 

In the case of a non-Abelian gauge symmetry with 
respect to which all standard model fields are singlets. 
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FIG. 2: Electron-paraphoton interaction induced in a renor- 
malizable model. 

gauge invariance forces the operators describing inter- 
actions of the new gauge bosons to include two field 
strengths. As a result, the new gauge bosons can be 
produced only in groups of two or more, and most op- 
erators have mass dimension eight or higher. The only 
exceptions are the operators analogous to the last two 
terms in Eq. Q), which lead to invisible Higgs decays. 

The operators shown in Eq. ^ describe physics only at 
scales below ~ M. However, simple renormalizable mod- 
els generate these operators after integrating out some 
heavy states. A generic feature of these models is the 
presence of fields charged under U{l)p. The lightest par- 
ticle of this type is stable and may be a viable dark matter 
candidate p] provided it does not carry color or electric 
charge. It turns out though that even if that particle 
is a singlet under J7(l)y, its electric charge cannot be 
exactly zero. This is because a kinetic mixing of the 
U{l)p X U{1)y gauge bosons is induced at some loop 
level, and the SL{2, R) and 5*0(2) transformations that 
diagonalize the kinetic terms shift the electric charges of 
all fields carrying U{l)p charge A dark matter parti- 
cle with mass of order 1 TeV must have an electric charge 
less than 10~^, because otherwise it would have left an 
imprint on the cosmic microwave background [12]. 

In models where the dimension-six operators arise at 
one loop, there is always a field charged under both U{l)p 
and U{1)y- A loop with this field in the internal line 
would induce a gauge kinetic mixing which in turn would 
lead to an electric charge for the dark matter particle 
that may be larger than 10^*. This problem is avoided 
in some models that generate the dimension-six operators 



at two loops: the loop-induced kinetic mixing is negligible 
if all fields charged under U{l)p have zero hypercharge, 
and the renormalizable tree-level kinetic mixing is absent 
if one of the C/(l)'s is embedded in a non-Abelian group. 

As an example, consider a new scalar, e, with the same 
gauge charges as cr, a gauge singlet Dirac fermion iV, 
and a eejiN coupling. If in addition there is a scalar 
and a fermion ip charged only under U{l)p, and a (j)^N 
coupling, then en couples to the paraphoton at two loops, 
as shown in Fig. 2. A Higgs Yukawa coupling inserted on 
an en external line generates the last operator in Eq. Q • 
If the new particles have masses of order M, and the 
new gauge and Yukawa couplings are of order unity, then 
Ce, c^, Cefj. are given by a two- loop factor of order (47r)~'*. 
Consequently, the limit on M is rather loose: M ^ 100 
GeV from Eq. (|14l) . The only [/(l)p-charged particles are 
(j) and V': and the lightest of them is a cold dark matter 
candidate. Further studies are necessary to determine 
the region of parameter space where the dark matter halo 
does not collapse too fast due to 7' emission. 

The quark-paraphoton moments are induced in this 
model only at three loops and arc negligible. Similar 
models, with e replaced by a scalar having the same 
charges as one of the quark fields, induce at two loops 
only the quark-paraphoton moments, and the tightest 
limit on AI in that case is set by 7' emission from su- 
pernova 1987A: M Z 40 GeV for c„ - i'in)'^- In the 
presence of an H^H(j)^(f> coupling, the last two operators 
in Eq. ^ are induced at one loop, and the gauge kinetic 
mixing arises at two loops. 

It is intriguing that massless gauge bosons other than 
the photon may interact with ordinary matter. The 
rather weak bound, below the electroweak scale in per- 
turbative models, on the scale M that suppresses such 
interactions makes it possible to search in collider ex- 
periments for the underlying dynamics that generate the 
dimension-six operators. It would also be interesting to 
investigate alternative experimental methods of search- 
ing for massless gauge bosons. 
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